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We present analytic expressions of the Coulomb interaction between a spherical and a deformed 
nuclei which are valid for all separation distance. We demonstrate their significant deviations from 
commonly used formulae in the region inside the Coulomb radius, and show that they remove various 
shortcomings of the conventional formulae. 
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I. INTRODUCTION 

The Coulomb interaction between two extended objects is a fundamental quantity in many problems of physics. One 
example is a heavy-ion collision, where the standard procedure is to approximate the Coulomb interaction between 
the projectile and target nuclei by simply replacing the target radius by the sum of the projectile and target radii or 
by the so called Coulomb radius in the formula for the Coulomb interaction between a point charge and an extended 
target nucleus Qj. This approximation works well as long as the physically relevant region is outside the Coulomb 
radius. However, one definitely needs to improve the formulae if the region inside the Coulomb radius becomes to 
play an important role. 

In this paper, we present analytic expressions of the Coulomb interaction between a spherical projectile and a 
deformed target which are valid for any separation distance between them and remove various shortcomings in the 
standard formulae. Since the formulae would have general values, here we mainly concentrate on derivation of the 
analytic expressions and comparing with commonly used formulae. The application to actual problems will be reported 
in separate papers. 

In sect. 2 we explain the basic idea of the method. It gives various components of the Coulomb interaction such 
as the bare Coulomb interaction or inelastic scattering form factors in the linear or higher order couplings in terms 
of a one-dimensional Fourier integral. In sect. 3, we present their analytic expressions obtained by computer assisted 
performance of the Fourier integrals. In sect. 4 we compare the bare potential and the linear as well as the second 
order coupling form factors calculated by our formulae with those by commonly used formulae by taking the reaction 
of 16 with 238 TJ as an example. We summarize the paper in sect. 5, where we briefly mention some possible systems 
where the improved formulae will become important. 

II. COULOMB INTERACTION IN THE FOURIER TRANSFORM REPRESENTATION 

Denoting the densities of the projectile and target by pp and px, the Coulomb interaction between them is given 



by 



V(R,a) = I dfi I d? 2 -=, - -pp{rV)p T {ri) (1) 

\R+f 2 -fi\ 



where R is the position vector of the center of mass of the projectile measured from that of the target nucleus and 
describes their relative motion, a represents the ensemble of intrinsic coordinates, which are implicit in pp and px- 
They are the deformation parameters in the collective model which we adopt in the following. The key idea is to 
express the same quantity by the following Fourier transform representation [Mpl , 
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For simplicity we assume a uniformly charged object with a sharp surface for both projectile and target, i.e. 

Pl (r) = p^Q(R l (n)-r) (3) 

where Q(x) is a step function and the index i refers to either the projectile or the target. The angle dependent radius 
is given in terms of the deformation parameters as 



The normalization condition then gives 



,(o) 



3Z t e 



Ni(a) 



with 



Ni(a) = 



(4) 



(5) 



(6) 



One usually chooses aoo to conserve the volume. In that case, Ni is one. In the following, we consider a spherical 
projectile and a deformed target. Accordingly, we remove the index i to distinguish the projectile or target from the 
deformation parameters and denote Rp^ and simply by Rp and Rt- The integration over in eq.(|^) can be 
easily performed, and one obtains 



V(R, a) = ^pfRl V Y^(Q R ) / dkj x (kR) 



Jl(kR P ) (T) 

kR P ^ K ' 



where 



M Z ] ^ = J dr JX (kr)Y^(n r )p T (r) 

(o) „ r kR(n r ) 



x 2 j x (x)dx 



(7) 

(8) 
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By expanding with respect to the deformation parameters, one obtains 

(0) „x 

M £ ) ( fe ) = ^-[ SxoV^J x 2 JO (x)dx 



+ x 3 jx(x)axfj, 
+ x 2 {xj x (x) + 



x 2 dj\(x ) 
2 dx 



(10) 



where x = kRx- The first, second and third terms in the square brackets of eq.(|l0|) give the bare Coulomb interaction, 
the linear and the second order Coulomb couplings, respectively. We thus represent the total Coulomb interaction as 



V(R, a) = V {0) (R) + V (1) (R, a) + V {2) (R, a) + 

where the first three terms are given by 

V (0) (R) = Z P Z T e 2 F {0 \R) 

vW(R,a) = Z P Z T e 2 Y J F { x \R)Y^{n B )a Xli 

V^(R,a) = Z P Z T e 2 £ £ F® (R)Y^(n R )a Xl ^a X24l2 

Ai ,A2 ifJ-2 A,/x 
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with the form factors defined by 
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(15) 
(16) 
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If we take the rotating frame where the z-axis is chosen to be parallel to the coordinate of the relative motion R, as is 
often done in the studies of heavy ion fusion reactions, and if we assume an axially symmetric quadrupole deformation 
for the target nucleus, then the angular momentum algebra can be explicitly carried out and we obtain the following 
expressions for the linear and the leading second order Coulomb couplings 



vW(R,fo,0) = Z P Z T e 2 Y, fPW\Y*d(9,0) 



(18) 



A=2,4,6 



V^(R,(3 2 ,9) = Z P Z T e 2 \Fl%(R)^-LY w (0,0) + FQ(R)~^Y iO (0,O) 
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(19) 



where 02 is the quadrupole deformation parameter of the target nucleus, and 9 the Euler angle to specify the orientation 
of its axially symmetric axis in the rotating frame. 

The important achievement of the Fourier transform method is that one needs to perform only one dimensional 
integral, whose results we will present in the next section. Before moving, we wish to comment that one can easily 
extend the same procedure to incorporate the surface diffuscness of the colliding nuclei by smearing the sharp surface 
in eq. (^) with a Yukawa function . The change is simply that the integrand of the Fourier transform representation 
of each form factor eqs.(|l5|) through (|l7]) gets two additional factors corresponding to the Fourier transforms of the 
Yukawa function which specify the surface properties of the projectile and target nuclei. 



III. ANALYTIC EXPRESSIONS OF THE BARE COULOMB POTENTIAL AND COUPLING FORM 

FACTORS 



The Fourier integrals on the r.h.s. of eqs. @ through (0) are tedious but still doable analytically by hands, since the 
integrands are given only by polynomials and trigonometric functions. However, computer assisted analytic integration 
is much easier and practical. Introducing the parameters Rc and Rcc by Rc = Rt + Rp and Rcc = \{Rt — Rp)\i 
the results read, 



(1) bare Coulomb interaction, 

r vr, 



V (0) (R) = Z P Z T e { 



{—R 5 + 24(R P + R T ) 3 (-R 2 P + 3R P R T - R 2 T ) + A5R(-R 2 P + R 2 ,) 2 
+15R 3 (R 2 P + R 2 T ) - WR 2 (Rp + i^)}/(160i?|i4) 
+(R% - m%Rl + im%R% - m 2 P R% + R b T )/(32RR 3 P R 3 T ] 



I {-5R 2 + 3{-R 2 P + 5R?t)}/(10R 3 r ), 

These formulae agree with those derived in Q and used in ||. 
(2) quadrupole Coulomb coupling form factor of linear order, 
3i?§,/(5R 3 ), 



(R > Rc) 

(Rc> R> Rcc) 
(Rcc > R > 0) 



3R 2 /(WR T ) + {3R 5 - \2R 3 (3R 2 p + R 2 ,) + im(-3R% 
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+3(R% 
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R%)}/{2bm 3 p R 3 T ) 
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(R > Rc) 

(R c >R> Rcc) 
(Rcc > R > 0) 
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(3) hexadecapole Coulomb coupling form factor of linear order, 
R T /(3R 5 ), 
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(R) 



(R > Rc) 



18R 5 (7R 2 P 



18R 
10R e D R: 



Pip 



R 4 /(6R T ) + {7R 7 
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+ 9R 3 (-35R P + WR 2 P R T + Rip) 
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(4) hexacontatetrapole Coulomb coupling form factor of linear order, 

( 3R 6 T /(13R 7 ), 



3fl 6 /(26fl^) + {99R° - 216R 7 (11R 2 P + R 2 T ) + 84R 5 (-99R 4 P + 18R 2 P R 2 T + R%) 



(R c > R> Rcc) 
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+24R 3 (231R e P - 189R 4 P R 2 T + 21R 2 P R 4 T + R b T ) + 18R(-231R% + 429R P R 2 T 
-2WR 4 pR T + 29R 2 P R ( T + R% )}/ (65536 RpR 7 T ) 

+3(99flJ? - 315K^K^, + 350R e P R^ - 150R 4 P R%, + 15R 2 P R S T + R]?)/ (8192 RR 3 pR 7 T ) 
+21(-URp + 54RpRp - W5R 8 P Rp + WORpRp - 45R 4 P R% + 6R 2 P Rp + R%? )/ (16384R 3 R 3 r Rp) 
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(23) 



(Rc > R > Rcc) 



(Rcc > R > 0) 



(5) quadrupole Coulomb coupling form factor of second order, 

' 6R 2 T /(5R 3 ), 



r (2) 
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(R > Rc) 



3R 2 /(20R 3 „) - {3R 5 - 12R 3 (3R 2 P - R 2 T ) - 18R(3R P + 2R 2 P R 2 T + 3R%, )} / (512RpR%) 
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+3(3R% + 29R 6 p R? r - 270R%Rp + 512R 3 P R% - 300R 2 P R% + 35 R%) / (25m R 3 R 3 P R%) 



-3R 2 /(10R%.), 
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(Rc > R> Rcc) 
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(6) hexadecapole Coulomb coupling form factor of second order, 
( R 4 T /(R 5 ), 
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IV. COMPARISON WITH COMMONLY USED FORMULAE 



We now compare our formulae with three commonly used models, which are given by 
(A) model I (point projectile model) 



4 



vi o) {R)= ZpZr£. (26) 
it 

3 Rn 



F "'<«> = 2ATT Jir <»> 



(B) model II (uniform charge model 1) 



(29) 



V(V(R) - Z P Z T e 2 | ;^'i?c)[3 - (^) 2 ], (* J ifc) 

[R) -^WTT)\R x /{Rc) x+ \ (R<Rc) m 



with 



(C) model III (uniform charge model 2) 



V<°>(J2) - Z P Z T e 2 | 1/^) [3 _ ( _g_ )2]) ( * J ^ (32) 

t m - 3 / Rt/R x+1 , (R>Rt) 

(H) - (2A + 1) \ R x /(R T ) x+ \ (R<Rt) 



The model I is used in almost all analyses of heavy-ion fusion reactions at energies near and below the Coulomb 
barrier. The model II can be often found in textbooks of heavy-ion collisions Q. Model I and II can be too crude 
concerning the bare potential. The model II and III have a shortcoming that each coupling form factor has a cusp at 
the separation distance, where the asymptotic formulae are matched to the formulae in the short distance region. The 
derivative of the coupling form factor is discontinuous at that distance. In writing eqs.(|30|) and (|3l|), we have assumed 
that the target nucleus has an axially symmetric static deformation with the Coulomb deformation parameter 
Eq.(|3l|) is the scaling condition of the deformation parameter to guarantee the correct coupling in the asymptotic 
region. 

Figs.l through 3 compare the bare Coulomb potential and the linear and quadratic coupling form factors F^} 2 4 6 (R) 
and i^_ 2 i(R) calculated by these formulae and by our improved formulae. The form factors have been multiplied 
with ZpZre 2 to make the ordinate of all figures have the dimension of energy. They have been calculated for the 
scattering of 16 O with 238 U. The Coulomb radius parameter has been chosen to be 1.06 fm. These figures clearly 
show the shortcomings of all the three commonly used simple models, which are solved by our new formulae. 

The important question is whether the deviations of the conventional models from our improved formulae play some 
significant roles in actual problems of physics, e.g. in analyzing heavy-ion fusion reactions at energies near and below 
the Coulomb barrier. In order to have some idea on this question, we again consider the fusion reaction of a spherical 
projectile 16 with an axially symmetric deformed target 238 U. It is now well known that the excitation of the ground 
state rotational band of 238 TJ plays an important role in enhancing the fusion cross section in this reaction. Instead of 
performing full coupled-channels calculations to take this effect into account, one often describes this reaction based 
on the no-Coriolis and the sudden tunneling, i.e. degenerate spectrum limit, approximations. In this case, the fusion 
probability for each partial wave J is given by first calculating the fusion probability Pj (E, 6) for a fixed orientation 
of the target nucleus 9 and then taking average over 6 

Pj(E) = - [ Pj(E,6)sm6d6, (34) 
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The effective potential for each orientation 9 is given by 

h 2 



Vj(R,I3 2> 9) = V N (R 7 f3 2 ,9) + ^-^J(J + l) + Z P Z T e 2 (F^(R) + (3 2 F { x ]} 2 (R)Y 2 o(9 7 0)) (35) 
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where Vn(R, P2, 0) is the nuclear potential and /x the reduced mass between the projectile and target. Since our 
aim is not to perform quantitative analyses, but to illustrate under what circumstances our improved formulae show 
its power, we took a simple model which considers only quadrupole deformation — 0.289 || for 238 U and linear 
coupling. The Gauss integral in eq.(|34|) is replaced by the Jn t?" +1 points Gauss quadrature if the rotational excitation 
is truncated at J m ax 

i 

where Qi are the angles, where the Legendre polynomial Pj max +i{0) becomes zero. 

The active angles are the zeros of P^cosO), which are about 30.55 and 70.12 degrees, if J max =2. In Fig. 4 we 
show the effective s-wave potentials for these two angles calculated by our improved formulae (the solid line) and 
by three conventional models (the dotted, dashed and dot-dashed lines). We assumed a Wood-Saxon potential 
for the nuclear potential, where the radius parameter has been chosen such that Rn = RjP + Rt^2Y2o(9, 0) with 
R$ = R P + R T = 1.06(Ap /3 + A}/ 3 ). 

We see clear differences among four models. However, these differences do not have any physical significance for 
16 + 238 U fusion reactions at low energies, since the deviation is localized well inside the barrier region, while the 
fusion probability is governed by the barrier property, which is the same for all four calculations. Similar situation 
will hold in general for medium weight heavy-ion collisions. This is a natural consequence of taking the same values 
for the nuclear and Coulomb radii. We can, however, think of some cases where the difference plays an important 
role as will be mentioned below. 



V. SUMMARY AND DISCUSSIONS 



We have derived analytic expressions of the Coulomb interaction, which are valid for any separation distance 
between the projectile and target, and have demonstrated their significant deviations from commonly used models. 
Our new formulae solve the cusp and discontinuity problems in the form factors and in their derivatives in commonly 
used models. 

We argued that these deviations will not cause any significant effects on the fusion reactions between two heavy 
nuclei such as the 16 + 238 U fusion reactions at energies near and below the Coulomb barrier, which have been very 
popular subjects of nuclear physics in the past decades. However, we have to keep in mind that this conclusion has 
been drawn by assuming the same values for the nuclear and Coulomb radii. It is then natural that the deviation of 
the conventional models from our improved formulae takes place well inside the barrier region. One interesting system 
will therefore be the system where the Coulomb radius is larger than the nuclear radius. One such example could 
be heavy-ion collisions induced by unstable neutron deficient isotopes. Another interesting case will be the fusion as 
well as elastic and inelastic scatterings between light heavy ions, where the absorption in the internal region is not so 
strong. In this case, the differences among 4 models at short distances will lead to quite different cross sections from 
each other. 
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Figure Captions 



Fig.l 

Comparison of bare Coulomb interaction. The solid, dotted, dashed and dot-dashed lines have been calculated based 

on our improved formulae, Model I through III, respectively. 

I ; ig-2 

Comparison of the linear Coulomb coupling form factor plotted in energy scale (see text). Notation is the same as in 
Fig.l. (a) quadrupole coupling (b) hexadecapole coupling (c) hexacontatetrapole coupling. 

Fi s- 3 

Comparison of the second order Coulomb coupling form factor plotted in energy scale. The solid and the dotted lines 
are for our improved formulae and Model I. 

Fig.4 

Comparison of the effective s-wave barrier calculated based on eq. (|3^) with our improved formulae and 3 conventional 
models. The results are shown for two orientations, 30.6 and 70.1 degrees. The lower and the upper barriers correspond 
to the former and the latter, respectively. 
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Fig. 1 
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